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We report on the progress of our ongoing project to quantify the role of the charm quark in the 
non-leptonic decay of a kaon into two pions. The effect of its associated mass scale in the dy- 
namics underlying the AI = 1/2 rule can be studied by monitoring the dependence of kaon decay 
amplitudes on the charm quark mass using an effective AS = 1 weak Hamiltonian. In contrast to 
commonly used approaches the charm quark is kept as an active degree of freedom. Quenched 
results in the GIM limit have shown that a significant part of the AI = 1/2 enhancement is purely 
due to low-energy QCD effects. Moving away from the GIM limit involves the computation of 
diagrams containing closed quark loops which requires new variance reduction techniques in or- 
der to determine the relevant weak effective low-energy couplings. We employ a combination 
of low-mode averaging and stochastic volume sources in order to compute these diagrams and 
observe a significant improvement in the statistical signal. 
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1. Introduction 

The decay of a neutral kaon into two pions in a state with isospin / is described by the transition 
amplitudes 

iAje i5 ' = ({nn)!\H w \K°) , (1.1) 

where H w is the AS = 1 effective weak Hamiltonian and 5/ the 7T7r-scattering phase shift. From 
experiments it is deduced that the kaon decay amplitude into two pions with total isospin / = is 
about twenty times larger than the amplitude into a state with 1 = 2, i.e. 

The observed enhancement of Aq, referred to as AI = 1/2 rule, remains one of the long-standing 
problems in hadron physics. Within the Standard Model, short-distance QCD and electroweak ef- 
fects yield only a moderate enhancement. Therefore, the main contribution is expected to come 
from long-distance, i.e. non-perturbative, QCD effects or, if this is not the case, from new physics. 
Lattice QCD which is the non-perturbative approach to the gauge theory of the strong interaction 
through regularized, Euclidean functional integrals is the only known technique that allows to at- 
tack the problem from first principles, and possibly to reveal the origin of the AI = 1/2 rule. 
In the low-energy regime of QCD various sources for the enhancement are possible. These include 
pionic final state interactions at around 100 MeV; physics at an intrinsic QCD scale of Aqcd ~ 250 
MeV; or physics at the scale of the charm quark, i.e. around 1.3 GeV. It remains unclear whether 
the experimental observation is the result of an accumulation of several effects, or mainly due to a 
single cause or mechanism. 

A theoretically well-defined strategy to disentangle non-perturbative QCD contributions from var- 
ious sources was proposed in Ref. [1], with the specific aim to reveal the role of the charm quark 
in the explanation of the AI = 1/2 rule. The possibility that the enhancement is mainly due to its 
mass being decoupled from the light quark mass scale was pointed out a long time ago [2]. 
In order to study non-leptonic weak decays of light hadrons one naturally turns to an effective 
theory given that the length scale of the weak interaction is two orders of magnitude smaller than 
the typical size of the involved hadrons. By means of the Operator Product Expansion (OPE) the 
non-local exchange of a gauge boson is replaced by an effective, local point interaction. Consid- 
ering QCD corrections the resulting effective Hamiltonian is given in terms of a sum over local 
fermion operators weighted by the Wilson coefficients. Thereby long and short distance physics 
are separated. The Wilson coefficients which include all high-energy effects are known at 2-loop 
of Perturbation Theory whereas the computation of the low-energy contributions to the transition 
amplitudes, i.e. the matrix elements of the operators, requires non-perturbative techniques. To this 
end the corresponding counterparts of the operators have to be constructed in the regularized lattice 
theory and put between a kaon and two pions in the initial and final state respectively. Despite the 
fact that many challenges arise in the direct computation of K — > nn remarkable progress [3] has 
been made for the past few years. 

In an alternative approach which alleviates the technical/computational challenges of computing 4- 
point functions in a large volume one considers the K — > n and K — > vacuum transitions, which are 
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then related to the physical ones by means of Chiral Perturbation Theory (ChPT) [4]. This forms 
the conceptual framework of the strategy proposed in Ref. [1] which is sketched in the following. 

2. AI = 1/2 rule on the lattice 

The crucial part of the proposed strategy consists in keeping the charm quark as an active 
degree of freedom, such that the theory has a softly broken SU(4)l x SU (4) r chiral symmetry 
where the up, down and strange quark have a degenerated mass at all stages. The role of the 
charm quark in the dynamics underlying the AI = 1/2 rule for kaon decays can then be studied by 
monitoring the dependence of kaon decay amplitudes on the charm quark mass m c . This is done 
in two steps: i) Set m c equal to the light quark masses, i.e. m u = mj = m s = m c (GIM limit), ii) 
Increase the charm quark mass towards its physical value. 

After the Operator Product Expansion to lowest order the AS = 1 effective weak Hamiltonian H w 
is given by 

H w = t {kfQt+kTQT}, (2.1) 

myi w ;=i 

where V qq > are CKM-matrix elements and kf 2 the Wilson coefficients, which incorporate all the 
high-energy effects. The four-quark operators Qf are given by 

Qf = { {sy^P.u){uy^P.d) ± (sy^P.d^uy^P.u)} - (w -> c). (2.2) 

Under SU(4)i the operator Qf(Qj) transforms as an irreducible representation of dimension 
84(20), while both Qf are singlets under SU(4)r. In the case of a diagonal mass matrix the 
operators Qf are of the form 

Qf = (m 2 u - m 2 c ) {m d (sP+d) + m s (sP-d)} , (2.3) 

where P± = \ /2(1 ± 75). Even though Qf do not contribute to the physical matrix elements, they 
are allowed by the underlying symmetries as a part of the effective Hamiltonian and mix with Qf 
under renormalization if m c ^m u . 

At leading order in ChPT, the ratio of amplitudes IA0/A2I is related to a ratio of low-energy con- 
stants (LECs) gf via 

2 8i ) 

Here, the LECs gf are the couplings multiplying the counterparts of the four-quark operators Qf 
in the effective Hamiltonian of the low-energy theory [4] . Obviously, the hierarchy of the ampli- 
tudes is directly related to a hierarchy in the couplings. The LECs can be determined by computing 
suitable correlation functions of Qf and Qf in LQCD and matching them to the corresponding ex- 
pressions in ChPT. The matching itself can be performed either in the standard p-regime of ChPT, 
or in the £-regime. The advantage of the latter is that no new LECs appear at next-to-leading order, 
which a priori may allow for a better control of systematic uncertainties. 

While in the first step (m c = mu g h t ) of the proposed strategy SU (4)-LQCD is matched to SU (4)- 
ChPT, in the second step (m c — > m p hy S ical) one matches the lattice results to SU (3) -ChPT after 
having integrated out the charm quark. 
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The complicated renormalization and mixing patterns of four-fermion operators usually encoun- 
tered in lattice formulations can be avoided through the use of Ginsparg- Wilson fermions [5], i.e. 
by considering Dirac operators that satisfy {D, 75} = aDy$D. Our choice is the Neuberger-Dirac or 
overlap operator [6], i.e. 



D 



1 



A = 1 + s — aD\y , a 



\s\ < 1. 



(2.5) 



Here, Dw refers to the Wilson-Dirac operator and the tunable parameter s allows to improve the 
locality properties of D. Introducing the modified quark field ¥ = (1 — aD/2)m guarantees that 
the renormalization and mixing of Qf are like in the continuum theory and, in particular, that no 
mixings with lower-dimensional operators with enhanced divergences occur [7]. The combined 
use of a SU (4)-flavor symmetry and chiral fermions, thus, leaves one with logarithmic divergences 
only. 

The use of dynamical overlap fermions is computationally very expensive, and the first studies 
(e.g. [8-10]) have been carried out in the quenched approximation. Intrinsic QCD contributions to 
the enhancement can be isolated by determining gf in the theory with m u = mj = m s = m c . In this 





Figure 1: Diagrams to be computed: "8"-diagram (left) and "eye"-diagram (right). 



case only the figure "8"-diagram (left of Figure 1) has to be dealt with. A moderate enhancement 
is observed [10], namely 




(2.6) 



This is not large enough to explain the experimental ratio but is already significant and cannot be 
attributed to penguin diagrams 1 . 

In the next step IA0/A2I has to be monitored as m c departs from the mass-degenerate limit towards 
its physical value. Thereby, the specific contribution of the charm quark to the AI = 1/2 rule can 
be investigated in detail. However, as soon as m c 7^ m u a new kind of diagram, referred to as 
"eye"-diagram (right of Figure 1), emerges and spoils the statistical signal. Its correlation function 
consists of two terms, a "color-connected" (con) and a "color-disconnected" (/fo) term which read 

C c ™ - (Tr {y^S(z,x)y Q P^S(x,y)y P^S(y,z)y^PS u/c (z,z)}) (2.7) 
Cg - (Tr {y^PS^zyyoP-Six^yoP-S&z)} Tr {y^S u/c (z,z)}) . (2.8) 

The challenge consists in the closed quark loop S u / C (z,z) with either a charm quark (c) or an up 
quark (u) running through it. By means of conventional techniques of computing quark propagators 

'Note, that the number of eq. (2.6) incorporates a factor 2-3 provided by high-energy effects calculated by means 
of Perturbation Theory. 
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involving point sources, it is not possible to sample over the coordinate z and noise dominates over 
the statistical signal. More sophisticated "all-to-all" propagators have to be developed to compute 
these contractions. 



3. All-to-all propagators 

At low quark masses the numerical computation of correlation functions is hampered by huge 
statistical noise which can be traced back to space-time fluctuations in the lowest eigenfunctions of 
the Dirac -operator [11, 12]. In the spectral representation the quark propagator can be written as 

V j Aj + m 

where v denotes eigenmodes and X the corresponding eigenvalues. If the quark mass m is of the 
same order as the gap, AA = — At, between two consecutive eigenvalues, the low-lying spec- 
trum is discrete and the lowest modes have a big weight in the sum. Space-time fluctuations in 
these eigenfunctions can lead to large fluctuations in observables whenever peaks in the wave- 
function happen to be near the fixed-point of the point-to-all propagator. Low-mode averaging 
(LMA) [11,13] has been shown to be a suitable tool to reduce these fluctuations. It amounts to sep- 
arating the nio W lowest eigenmodes from the rest, i.e. to truncating the sum over all eigenmodes. 
Consequently, the quark propagator decomposes into a "low"(S') and "high"(S*) part 

S(x,y) = S l (x,y)+S h (x,y) = 1 "£ v ' ( *>® v ? (y) + S h (x,y). (3.2) 

V i Aj+m 

The latter lives in the orthogonal complement of the subspace spanned by the n\ ow lowest modes. 
Decomposing the quark propagator into a low-mode and high-mode propagator results in the split- 
ting of the correlation function as well. So, a common three -point function which consists of 4 
propagators is split into 5 terms with in total 16 distinct contributions. Schematically, 

Q^pt (-4111 i ^-jlllh | gllhh | £<lhhh _^_^-<hhhh 3) 

# diagrams = 1 4 6 4 1 

where the upper indices / and h denote the number of low-mode and high-mode propagators accord- 
ing to the decomposition of eq. (3.2). The statistical signal for the correlation function is enhanced 
by exploiting translational invariance in the terms with at least one low-mode propagator. Since the 
latter is known for all space-time locations these contributions can be sampled over many different 
source points. 

Furthermore, to improve the statistical signal at fixed computational cost it turned out that instead 
of merely increasing n\ ow it is advantageous [14] to keep n\ ow reasonably small and construct "ex- 
tended" all-to-all propagators, which allows to sample also over the high part of the propagator. 
Here, the mode itself is used as a source for an additional inversion. More precisely, taking the 
source to be the left-projected eigenmode (7oP_)v,- at the fixed timeslice t = tf, the solution vector 
S ext reads 2 

$i ( x )\)'o=tf 



z j ^S h (x,f,tf)^P-) Vi (f,tf), (3-4) 



Note, that a semicolon ";" is used to write the spatial and temporal part of a 4-vector separately, i.e. y = y;yo- 
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where an average over the spatial coordinate y is automatically performed. LMA with extended 
all-to-all or conventional point-to-all propagators for the "high" part S h is crucial to obtain a signal 
in the £ -regime even for the simple case of the figure "8"-diagram. 

The "eye"-diagram, however, which is required when decoupling the charm quark mass is still 
dominated by statistical noise. To reduce its variance, stochastic volume sources (SVS) [15] and 
dilution [16] techniques are used to estimate S h stochastically. That is, low-mode averaging is com- 
bined with stochastic sources (LMA + SVS). The resulting all-to-all propagators allow to average 
over the intersection point of the 4-quark lines, i.e. the point where two quark propagators attach 
to the closed quark loop. 

To compute stochastic all- to-all propagators an ensemble of r = l,...,N T random noise vectors, 
{r]^(xo,x)}, is generated for each gauge configuration. These source vectors are created by as- 
signing independent random numbers to all components, i.e. to all lattice sites, color and Dirac 
indices and have to obey the following two conditions 3 

(n a a (x ,x)) SK ee lim 1 £ (n^YaM = (3.5) 

Va(^o)(Tl%(y,yo))^ = 8 xm 5^8 a p8 ab . (3.6) 
Then one can invert for each of these noise vectors and obtain an estimate for the full propagator 

;«&W(i? t )J(y)) Bc = sS(^y), (3-7) 

where the individual solution for each of the N r noise vectors is given by 

(* w )iw=EE^.«)(»? (r) );w- (3-8) 

For the method to be efficient, one must ensure that the intrinsic stochastic noise does not over- 
whelm the gain in information provided by having access to the entire propagator matrix. An es- 
sential step towards a significant variance reduction is taken by restricting the support of the source 
vector to individual timeslices, Dirac or colour components [16]. Such "dilution schemes", and, in 
particular, the so-called time-dilution are widely used in the computation of hadronic properties. In 
the time-dilution scheme the non-zero components of the random source vector are restricted to a 
single timeslice, e.g. yo = 

n(y) = {]> i9) > if l° = °' • (3.9) 

,KJ |0, otherwise 

In this work dilution is applied in spin, color and time, i.e. each of the N r noise vectors of the 
ensemble has random entries only on one spin-color component of a single timeslice with all other 
entries set to zero. 

4. Results 

A single quenched lattice is used to study the performance of LMA + SVS where volume 
sources are used for estimating S h stochastically. The bare coupling constant is /3 =6/gl = 5.8485 

3 Latin(Greek) letters denote color(spin) components. 



6 



Contribution of the charm quark to the AI — 1/2 rule 



Eric Endress 



0.15 
0.1 
0.05 


-0.05 
-0.1 
-0.15 



I 



i i i i r 



» LMA + SVS 
» LMA 





0.2 



0.1 



-0.1 



light 



1 


1 1 




1 






I 






— I 




n 




\ 






1 


1 1 




1 



heavy 



£,aU ^hhhh ^lhlili ^llhh „Ulh £,1111 



£,all ^.hhhh lhhh ^Uhh lllh 1U1 



Figure 2: Individual contributions to the ratios R± classified by their number of low-mode propagators for 
a light (top) and heavy (bottom) charm quark. Blue data points (displaced to the right for better visibility) 
result from LMA only, whereas LMA + SVS is used for the red data points (the stochastic estimate is only 
considered for a diagram if its noise is reduced, otherwise the LMA data is kept). The contribution with 
low-mode propagators, i.e. C , is not computed stochastically and C" 11 refers to the sum of all 5 terms. 



which corresponds to a lattice spacing a ~ 0. 12 fm, and a volume Va~ 4 = 16 3 x 32. The light bare 
quark mass is am;» = 0.02, resulting in a pion mass m n rj 320 MeV. Two charm quark masses 
are considered: am c = 0.04 = 2 x amug^ and am c = 0.2 = 10 x araj,^. Twenty low-modes are 
computed for each of the 200 quenched configurations. The volume sources are diluted in time, 
spin and color. Stochastic estimators are used for the loop propagator and, if required, also for the 
other propagators in order to be able to average over the position of the 4-quark line intersection. 
By applying LMA the "eye"-diagram splits into 5 distinct terms grouped by the number of low- 
mode propagators according to eq. (3.3). In the following ratios of these terms are shown. More 
precisely, the ratios are defined by 



fl=F(|*0-Zo|,bo-Zo|) 



C%{\xq - zol bo - zoQ ±gg (1*0 - zo\, bo - ZoQ 
C 2 (*o)C 2 (3'o) 



(4.1) 



where C 2 is the two-point function of the left-handed current 7o = fi'YoP-^)- 
The overall improvement of LMA+S VS compared to LMA is illustrated in Figure 2. It reveals that 
LMA+SVS is effective for the terms consisting of 2 or 3 low-mode propagators. The variance is 
reduced significantly, most notably when the charm quark in the closed loop is heavy. In the latter 
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case the absolute error of the sum of the 5 terms is roughly halved. For the term with a single low- 
mode propagator the technique of LMA+SVS shows no significant improvement; presumably in 
this case the use of multiple independent stochastic estimates for several high-mode parts increases 
the intrinsic stochastic noise and the technique deteriorates its performance. 

5. Summary and outlook 

We have reported on the progress of our ongoing project to understand the role of the charm 
quark and its associated mass scale in non-leptonic decays of kaons into two pions. When the 
charm quark mass is decoupled from the light quark masses, it is hard to obtain statistical signals for 
"eye"-diagrams. A combination of low-mode averaging and stochastic volume sources is applied 
to cure this. We observe a significant variance reduction for several contributions, even though the 
overall error remains sizable. In the next step of this project the results for the bare ratios will be 
renormalized. To this purpose the contributions of the operators Qf have to be taken into account. 
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